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Abstract 

The correlator of topological charge densities x(Q 2 ) m QCD is calculated in 
the domain of low Q 2 . Basing on the Ward identities the low energy theorems 
are proved. The value of the first derivative x'(0) found recently by QCD sum 
rules is used. The contributions of pseudoscalar quasi-Goldstone bosons as 
intermediate states in the correlator are calculated with the account of 1] — tt 
mixing. 
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I. INTRODUCTION 

The existence of topological quantum number is a very specific feature of non-abelian 
quantum field theories and, particularly, QCD. Therefore, the study of properties of the 
topological charge density operator in QCD 

cy ~ 

Q*W = ^ G U^)GUx) (1) 
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and of the corresponding vacuum correlator 

X (q 2 ) = i J d'xe^(0 | T {Q 5 (x), Q 5 (0)} | 0) (2) 

is of a great theoretical interest. (Here G™ u is gluonic field strength tensor, G^ u = 
(l/2)e flv x a G\ a is its dual, n are the colour indeces, n = 1, 2, ...N% — 1, N c is the number of 
colours, N c = 3 in QCD). The existence of topological quantum numbers in non-abelian field 
theories was first discovered by Belavin et al. [[J, their connection with non- conservation of 
£7(1) chirality was established by t'Hooft ||. Crewther Q derived Ward identities related 
to x(0); which allowed him to prove the theorem, that x(0) = in any theory where it is at 
least one massless quark. An important step in the investigation of the properties of x{q 2 ) 
was achieved by Veneziano || and Di Vecchia and Veneziano ||. These authors considered 
the limit N c — > oo. Assuming that in the theory there are Nf light quarks with the masses 
mi <C M, where M is the characteristic scale of strong interaction, Di Vecchia and Veneziano 
found that 

f Nf i r 1 

x(o) = (o|m|o> £ — , (3) 

V i m iJ 

where (0 | qq | 0) is the common value of quark condensate for all light quarks and the terms 
of the order rrii/M are neglected. [] The concept of #-term in the Lagrangian was succesfully 
exploited in in deriving of (3). Using the same concept and studying the properties of 
the Dirac operator Leutwyler and Smilga succeeded in proving eq.3 at any N c for the 
case of two light quarks, u and d. 

In this paper I calculate x{q 2 )i a t l° w I q 2 \^ M 2 . In ref.l the value (more precise, 
its nonperturbative part) was found. On the basis of QCD sum rules in the external fields 
the connection of with the part of the proton spin, carried by u, d, s quarks, S was 
established. By the use of experimental data on E, as well as from the requirement of 
selfconsistency of the sum rule, it was obtained: 

X'(0) = (2.3 ± 0.6) x 1(T 3 GeV 2 (4) 



The definition of x(q 2 ) used above eq.2, differs by sign from the definition used in Q- Q. 
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in the limit of massless u, d and s quarks. At low q 2 the terms, proportional to quark masses, 
are related to the contributions of light pseudoscalar mesons as intermediate states in the 
correlator (2). These contributions are calculated below. In such calculation for the case of 
three quarks the mixing of 7r° and r\ is of importance and it is accounted. 

The presentation of the material in the paper is the following. In Sec. II the low energy 
theorems related to x(0) are rederived with the account of possible anomalous equal-time 
commutator terms. (In [|J- || it was implicitly assumed that these terms are zero). In 
Sec. Ill the case of one and two light quarks are considered. It is proved, that the mentioned 
above commutator terms are zero indeed and for the case of two quarks eq.3 is reproduced 
without using iV c — > oo limit and the concepts of #-terms. In Sec. IV the case of three u, d, s 
light quarks is considered in the approximation m u ,md <C m s . The problem of mixing of 
7T° and 7] states |||§ is formulated and corresponding formulae are presented. The account 
of 7r — i] mixing allows one to get eq.(3) from low energy theorems, formulated in Sec. II for 
the case of three light quarks. (At m u ,md <C m s it coincides with the two quark case). In 
Sec.V the g 2 -dependence of x{q 2 ) was found at low | q 2 | in the leading nonvanishing order 
in q 2 /M 2 as well as in m q /M. 

II. LOW ENERGY THEOREMS. 

Consider QCD with Nf light quarks, m,; < M ~ 1 GeV, i — 1, ...Nf. Define the singlet 
(in flavour) axial current by 




(5) 



and the polarization operator 




(6) 



The general form of the polarization operator is: 




(7) 



Because of anomaly the singlet axial current is nonconserving: 




(8) 
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where Q$(x) is given by (1) and 

D(x) = 2i^2m i q i (x)'y 5 qi(x) (9) 

i 

It is well known, that even if some light quarks are massless, the corresponding Goldstone 
bosons, arising from spontaneous violation of chiral symmetry do not contribute to singlet 
axial channel (it is the solution of U(l) problem), i.e. to polarization operator P^ u (q). P^q 2 ) 
also have no kinematical singularities at q 2 = 0. Therefore 

P,u(q)q,qu = -P L (q 2 )q 2 (10) 

vanishes in the limit q 2 — > 0. Calculate the left-hand side (lhs) of (10) in the standard way 
- put q^q u inside the integral in (6) and integrate by parts. (For this it is convenient to 
represent the polarization operator in the coordinate space as a function of two coordinates 
x and y.) Going to the limit q 2 — > we have 

fim P^{q)q»q» = * I d 4 x(0 | T{2N f Q 5 {x), 2N f Q 5 (0) 

q 2 ^0 J 

+2N f Q 5 {x), D{0)+D(x), 2N f Q 5 (0)+D(x), D(0)} | 0) 

N f 

+4^771^(1 gi(0)g<(0) | 0) 

i 

+ / d 4 x(0 | \j 05 (x), 2N f Q 5 (0) ] | 0)S(x ) = (11) 

In the calculation of (11) the anomaly condition (8) was used. The terms, proportional 
to quark condensates arise from equal time commutator [jos^x), D(0)] xo=0 , calculated by 
standard commutation relations. Relation (11) up to the last term was first obtained by 
Crewther M. The last term, equal to zero according to standard commutation relations and 
omitted in J3|- 0, is keeped. The reason is, that we deal with very subtle situation, related 
to anomaly, where nonstandard Schwinger terms in commutation relations may appear. (It 
can be shown, that, in general the only Schwinger term in this problem is given by the last 
term in the lhs of (11): no others can arise.) Consider also the correlator: 

P,(q) =ijd i x e^(0 | T{j^(x), Q 5 (0)} | 0) (12) 
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and the product P tl (q)q ll in the limit q 2 — > (or q 2 of order of the m 2 , where m n is the mass 
of Goldstone boson). The general form of P^{q) is P^{q) = Aq^. Therefore nonvanishing 
values of P^q^ in the limit g 2 — ^ (or of order of quark mass m, if | q 2 |~ m 2 - this limit 
will be also intersting for us later) can arise only from Goldstone bosons intermediate states 
in (12). Let us us estimate the corresponding matrix elements 

(0 | | tt) = Fq, (13) 

(0 | Q 5 | tt> = F' (14) 

F is of order of m, since in the limit of massless quarks Goldstone bosons are coupled only 
to nonsinglet axial current. F' is of order of m 2 /^ ~ m, where f n is the pion decay constant 
(not considered to be small), since in massless quark limit, the Goldstone boson is decoupled 
from Q§ . These estimations give 

q 2 

P»% 2 m2 ( 15 ) 

q z — 

and it is zero at q 2 — > 0, m 2 7^ and of order of m 2 at q 2 ~ m 2 . In what follows I will 
restrict myself by the terms linear in quark masses. So, I can put P ll (q)q ll = at q — > 0. 
The integration by parts, in the right-hand side (rhs) of (12) gives: 

lim PM)% = - / d 4 x(0 I T{2N f Q 5 {x), Q 5 (0) + D(x), Q 5 (0)} \ 0) 

- J d 4 x(0 I [ j 05 (x), Q 5 (0) ] I 0)6(x ) = (16) 
After the substitution of (16) in (11) arise the low energy theorem: 

i J d 4 x(0 I T{2N f Q s (x), 2N f Q 5 (0)} \ 0) 

N f 

-1 J d 4 x{0 I T{D(x), D(<d)} I 0) -4j2 mi (0 I g,(0)g,(0) | 0) 

i 

+ i J d*x(0 I [ f 05 (x), 2N f Q 5 (0) ] I 0)S(x ) = (17) 
The low energy theorem (17), with the last term in the lhs omitted, was found by Crewther 



III. ONE AND TWO LIGHT QUARKS. 



Consider first the case of one massless quark, Nj — 1, m — 0. This case can easily be 
treated by introduction of #-term in the Lagrangian, 

The matrix element (0 | Q5 | n) between any hadronic state | n) and vacuum is proportional 

J d 4 x(0 I Q 5 {x) I n) ~ (0 I — InZ \ n) e=0 , (19) 

where Z = e lL and L is the Lagrangian. The gauge transformation of the quark field 
ip' — > e iai5 il) results to appearance of the term 

6L = a [ - (a s /4-K)G^G^ ] (20) 

in the Lagrangian. By the choice a = 9 the (9-term (18) will be killed and (d/d9)lnZ = 0. 
Therefore, x(0) = (Crewther theorem). The first term in (17) vanishes, as well as the 
second and third, since m = 0. From (17) we have, that indeed the anomalous commutator 
vanishes 

<0| [j 05 (x), Q 5 (0) ] X0=Q I 0) = 0, (21) 

supporting the assumptions done in 0]- [|]. 

Let us turn now to the case of two light quarks, u, d, Nf = 2. This is the case of real 
QCD, where the strange quark is considered as a heavy. Define the isovector axial current 

4? = (^7sM - dl^ 5 d)/V2 (22) 

and its matrix element between the states of pion and vacuum 

(0 I $ I 7T> = Uq^ (23) 

where g M is pion 4-momentum, f n = 133 MeV. Multiply (23) by q^. Using Dirac equations 
for quark fields, we have 
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2i - i 

-^=(0 | m u uj 5 u - m d d^ 5 d | vr> = — ^=<0 | (m u + m d )(uj 5 u - d~f 5 d) 

+ (m u + m d )((u^u - d~f 5 d) \ ir) = f v m%, (24) 

where m u , rrid are u and d quark masses. The ratio of the matrix elements in lhs of (24) is 
of order 

" (25) 



(o 


wy 5 u + d^ 5 d 


7T> 


(o 


— d'jsd 


tt) 



M ' 

since the matrix element in the numerator violates isospin and this violaton (in the absence of 
elecntomagnetism, which is assumed) can arise from the difference m u — m d only. Neglecting 
this matrix element we have from (24) 

4=(0 I u l5 u - d l5 d I tt) = J^L- (26) 
V 2 m u + m d 

Let us find x(0) from low energy sum rule (17) restricting ourself to the terms linear in quark 
masses. Since D(x) ~ m, the only intermediate state contributing to the matrix element 

J d 4 x(0 I T{D{x), D(0)} I 0) (27) 

in (17) is the one-pion state. Define 

D q = 2i(m u u-f 5 u + m d d-f 5 d). (28) 

Then 

(0 I D q I 7r) = i(0 I (m u + m d )(uj 5 u + dj 5 d) + (m u - m d )(uj 5 u - dj 5 d) \ n) 

= v^^M4, (29) 
m u + m d 

where the matrix element of singlet axial current was neglected and (26) was used. The 
substitution of (29) into (27) gives 
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d*xe^(0 | T{D q (x), D q (0)} I 0)^o = lim < ~— -j2( ) />* 

= -4 (m "- mrf)2 (0 | gg | 0) (30) 
m u + m d 

In the last equality in (30) Gell-Mann-Oakes-Renner relation [10] 

1 f 2 m 2 

(0 | qq | 0) = -- -J^- (31) 

2 m u + m d 

was substituted as well the £77(2) equalities 

(0 | uu | 0) = (0 | dd | 0) = (0 | qq | 0). (32) 

^From (17) and (30) we finally get: 

X (0) = i I d*x(0 | T{Q 5 (x), Q 5 (0)} | 0) = ^ (0 | gg | 0) (33) 
7 m u + m d 

in concidence with eq.3. 

In a similar way matrix element (0 | Q5 | n) can be found. Consider 

(0 I i M5 I tt> = Fg M (34) 



The estimation of F gives 



F iw - A (35) 



and after multiplying of (34) by q^ the rhs of (34) can be neglected. In the lhs we have 

(0 I D q I tt) + 2^,(0 I Q 5 I tt> = (36) 

The substitution of (29) into (36) results in 

(0 I Q 5 I tt) = --^^K (37) 
2a/2 m„ + m d 

The relation of this type (with a wrong numerical coefficient) was found in ||, the correct 
formula was presented in [|ll|]. From comparison of (33) and (37) it is clear, that it would be 
wrong to calculate x(0) by accounting only pions as intermediate states in the lhs of (33) - 
the constant terms, reflecting the necessity of subtraction terms in dispersion relation and 
represented by proportional to quark condensate terms in (17) are extremingly important. 
The cancellation of Goldstone bosons pole terms and these constant terms results in the 
Crewther theorem - the vanishing of x(0), when one of the quark masses, e.g. m u is going 
to zero. 
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IV. THREE LIGHT QUARKS. 



Let us dwell on the real QCD case of three light quarks, u, d and s. Since the ratios 
m u /m s , md/m s are small, less than 1/20, account them only in the leading order. When 
the u and d quark masses m u and rrid are not assumed to be equal, the quasi- Goldstone 
states 7T° and 77 are no more states of pure isospin 1 and correspondingly: in both of these 
states persist admixture of other isospin |8||| proporional to m u — m^. (The violation of 
isospin by electromagnetic interaction is small in the problem under investigation |J and 
can be neglected. The 77' — 77 mixing is also neglected.) In order to treat the problem it is 
convenient to introduce pure isospin 1 and pseudoscalar meson fields tps and ips in 577(3) 
octet and the corresponding states | P 3 ), | P 8 ) PIPl- Then in the 577(3) limit 

(0 | 3 { S I P 3> = Uq, (38) 

(0 | j$ | Ps) = Uq» (39) 

where 

3$ = («7m75« + dl^ld - 2s7 M 7 5 s)/ VQ (40) 

and is given by (22). The states | P3), | Ps) are not eigenstates of the Hamiltonian. In 
the free Hamiltonian 

H = ^rhl(pl + ifnfyl + (P 8 | P 3 )(f3<f8 + kinetic terms (41) 

the nondiagonal term ~ ^3^8 is present ( m^. and in (41) coincides with and up 
to terms quadratic in (P 8 | P3)). The nondiagonal term was calculated in |J on the basis of 
PCAC and current algebra (see also 0) 

(Ps I Ps) = ^=m;^i^ (42) 
V3 m u + m d 

The physical ir and 77 states arise after orthogonalization of the Hamiltonian (41) 

I tt) = cosO I P 3 ) - sin6 I P 8 ) 

I r?) = sinO I P 3 ) + cose I P 8 ) (43) 
10 



where the mixing angle 9 is given by (at small 9) @,[|: 

(Ps | P 3 ) _ (A I A) 1 ml m u -m d 

u = ~ 2 2 2 — = ~~ 2 — ; — v 44 ^ 

In terms of the fields ip^ and </2 8 PCAC relations take the form ||; 

d»j$ =Mmltps + (Ps\P3)(ps) (45) 

94?=A(m^ 8 + (P 8 |P 3 )^ 3 ) (46) 

Our goal now is to calculate the contribution of pseudoscalar octet states to the second term 
in the lhs of (17). It is convenient to use the full set of orthogonal states | P 3 ), | P 8 ) as the 
basis. Use the notation 

D = D q + D s D s = 2im s s l5Sj (47) 

where D q is given by (28). The matrix element 

(0 I D q I P 3 ) = ^^Urnl (48) 

can be found by the same argumentation, as that was used in the derivation of (29). In 
order to find (0 | D q | P 8 ) take the matrix element of eq.45 between vacuum and | Ps) 

(0 I I P 3 ) = U{P 8 I P 3 ) (49) 

The substitution in the lhs of (49) of the expression for d^jj^j through quark fields gives 

(0 I - d^d I P 8 ) = - m " ■ md (0 I + d^d I P 8 )- 

m u + m d 

. [2 2 m u -m d 

In a similar way take matrix element of eq. (46) between (0 | and | P 8 ) and substitute into 
it (50). We get 
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75 (0 



m u + m d - 



(m u - m d f 
m u + m d 

2 2 r I "' " ' " <' 



(«7 5 m + d7 5 cf) - 4m s s7 5 s | P 8 ) 



3 " \m u + m d/ 
As follows from SU (3) symmetry of strong interaction 

(0 | + dj 5 d I P 8 ) = -(0 I S7 5 s I P 8 ) 

up to terms of order m q /M, which are neglected. From (51), (52) we find: 



(0 I D s I P 8 > = -\\-Um% 



1 (m u -m d ) 2 
4m s (m u + m d ) 



1 + 



m u m d 



m s (m u + m d ) 



(51) 



(52) 



(53) 



2 rn u m d 



(0| Dj |P 8 )=4y-/.< K + mJ2 
in notation (28), (47). When deriving (53) the SU(3) relation 



2 4 2 m 
= -m 



A _ l m u + m d 

m u + m d \ 4 m s 



(54) 



(55) 



was used. In (53) the small terms ~ m u /m s , m d /m s are accounted, because they are 
multiplyed by large factor m 2 rj . In (54) small terms are disregarded. The matrix element 
(0 I D s I P 3 ) can be found from (46). We have 



-j={0\D q -2D s \P 3 )=f T (P 8 \P 3 ) 



(56) 



The substitution of (48) and (42) into (56) gives 



(0 I D s I P 3 > = 



(57) 



Equations (48), (53), (54) and (57) allow one to calculate the interesting for us correlator 

i J d 4 x(0 I T{D(x), D(0)} I 0) (58) 

when the sets | P 3 )(P 3 | and | P 8 )(P 8 | are taken as intermediate states. But | P 3 ), | P 8 ) 
are not the eigenstates of the Hamiltonian, they mix in accord with (41). Therefore the 
transitions (P 8 | P 3 ) arising from the mixing term in (41) must be also accounted. There 
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are two such terms. The one corresponds to the transition (0 | D s \ P$) (P 3 \ D q | 0) and its 
contribution to (58) is given by 

hm{-2(0 | D s | Ps)^(Ps I Ps)^(Ps I D q | 0)} = 

= 2«f^^V (59) 

The other corresponds to the transition between two D s operators, where (P3 | P3) enter as 
intermediate state. This contribution is equal to: 



lim j-(0 I D s I P 8 )^-^(P 8 I p 3 )_!_<p 8 I p 3 )__!_<p 8 I p s I 0) 
g 2 ^o [ q 2 — m 2 q 2 — m% q 2 — m 2 



(60) 



It is enough to account only matrix elements, with D s operators, since they are enhanced 
by large factor m 2 . All others are small in the ratio m^/m 2 . 
Collecting all together, we get: 



i I d*x(0 I T{D(x), D(0)} I 0) = flml 2 ^ 
J 1 \ m u 



- m d 



m u m d 



+ m di 



| 2 m s f l m u + m d 
m u + m d \ 4 m s 



1 - 



1 m u -m d 
2m s (m u + m d ) 

2m., 



1 - 2 



(m u + m d ) 2 

m u m d 
m s (m u + m d ) 



m u + m di 



2\m u + m d/ 



c2 2 



9 / m u - m d 



m u + m d 2\m u + m d 



(61) 



The first term in the figure bracket in (61) comes from (0 | D q | P3) 2 , the second - from 
(0 I D q I P 8 ) x (0 I P 8 I P s ), the third - from (0 | D s | P 8 ) 2 , the last two terms are from 
(59), (60). Adding to (61) the proportional to quark condensate term 



4(m u + m d + m s )(0 | qq | 0) 



(62) 



in (17), we finally get for 3 quarks at m u ,m d <C m s 



% f d 4 x(0 I T{2N f Q 5 (x), 2N f Q 5 (x)} | 0) = 36 m " m<i (0 | qq | 0) 



(63) 
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and 

X(0) = -^^(0 I qq | 0), (64) 

since in this case 4Nj = 36. Eq.64 coincides with eq.3, obtained || in iV c — > oo limit. This 
fact demonstrates, that N c — > oo limit is irrelevant for determination of x(0) (at least for 
the cases of two or three light quarks and at m ui m d <C m s ). %(0) for three light quarks at 
m u , ma Cms - eq.64 coincides with x(0) in the two light quark case, (see [7] and [33]) i.e in 
this problem, when m u , rrid <C m s , there is no difference if s-quark is considered as a heavy 
or light - it softly appears in the theory. 

Determine the matrix elements (0 | Qs | tj) and (0 | Q§ \ tt). Following JTSfl consider 



(0 | J M 5 I V) = Fq, (65) 

F is of order of f n (m s /M) and can be put to zero in our approxamation. By taking the 
divergence from (65), we have 

(0 | D s + 6Q 5 | V) = (66) 

The use of (53) gives (the tt — r] mixing as well terms of order m u /m s , md/m s may be 
neglected here): 

(0 I Q 5 I V) = l\ll^ m v (67) 



Relation (67) was found in By the same reasoning it is easy to prove that 

(0 | D q | P 3 ) + (0 | D s | P 3 ) + 6(0 | Q 5 | P 3 ) = (68) 

The first term in (68) is given by (48), the second one is zero according (57). For the last 
term we can write using (43) 

(0 | Q 5 | P 3 > = (0 | Q 5 I tt) + 6(0 | Q 5 | P 8 ) (69) 

Eq.'s (68), (69) give 

(0 I Q 5 I tt) = -^Uml^^ (70) 
2a/2 ^« + m d 

- the same formula as in the case of two light quarks. 

It is clear, that the presented above considerations can be generalized to the case, when 
u,d and s-quark masses are comparable. The calculation became more cumbersome, but 
nothing principially new arises on this case. 
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V. Q 2 DEPENDENCE OF X {Q 2 ) AT LOW Q 2 . 



Let us dwell on the calculation of the g 2 -dependence of x(<? 2 ) a t l° w | <7 2 | in QCD, re- 
stricting ourselves by the first order terms in the ratio q 2 /M 2 , where M is the characteristic 
hadronic scale, M 2 ~ 1 GeV 2 . By x(<? 2 ) I mean its nonperturbative part with perturba- 
tive contribution subtracted from the total x{q 2 ) defined by (2). The reason is that the 
perturbative part is strongly divergent, its contribution would be strongly dependent on 
regularization proceedure and, therefore, physically meaningless. In this domain of q 2 x{q 2 ) 
can be represented as 

X{q 2 ) = X(0) + X\0)q 2 + R{q 2 ) ~ R{0) (71) 

x(0) for the QCD case-three light quarks with m u , rrid <C m s - was determined in Sec. IV. 
x'(0) (its nonperturbative part) for massless quarks was found in [0] basing on connection 
of x'(0) with the part of proton spin carried by u, d, s quarks. Its numerical value is given 
by (4). What is left, is the contribution of light pseudoscalar quasi-Goldstone bosons R(q 2 ), 
which has nontrivial g 2 -dependence and must be accounted separately. R{q 2 ) vanishes for 
massless quarks and did not contribute to x'(0), calculated in |]. R(0) must be subtracted 
from R(q 2 ) since it was already accounted in x(0). R(q 2 ) can be written as 

R{q 2 ) = "(0 | Q 5 | vr) 2 ^^ - (0 | Q 5 | V) 2 ^-^ (72) 

qz _ m z gz _ m z 

The problem of rj — ir mixing is irrelevant in the difference R{q 2 ) — -R(O) in any domain 
| q 2 |~ m 2 and | q 2 |~ m 2 . The matrix elements entering (72) are given by (67), (70). 
Taking the difference R{q 2 ) — -R(O) and using the Euclidean variable Q 2 = —q 2 , we have 



x(Q 2 ) = x(0) - x'(0)Q 2 - l -flQ 2 



'm u -m d \ 2 ml 1 m 2 



(73) 



\m u + m d J Q 2 + ml 3Q 2 + m 2 _ 

Eq.(73) is our final result, where x(0) is given by (64) and x'(0) by (4). The accuracy of 
(73) is given by the parameters Q 2 /M 2 , m^/M 2 , m 2 /M 2 1, (two last characterize the 
accuracy of SU(3) x S77(3)). At Q 2 w m 2 the last term comprise about 20% of the second 
(the first term is very small, x(0) ~ —4 - 10 -5 GeV A ). Evidently, the last term is much bigger 
in the Minkovski domain, Q 2 < 0, since there are pion and rj poles. As was mentioned above, 
x(0) found here concides with x(0) obtained in |J by considering large N c limit. However, 
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the g 2 -dependence is completely different. Namely, for x(<? 2 ) i R @ was found the relation 
(eq.(A4') in §) 



x{q 2 



aFl 
~2N r 



AT 



N, 



where the Goldstone boson masses fi 2 are related to quark condensate by 

tf = -2m~(0\qq\0) 



(74) 



(75) 



and a is some constant of order of hadronic mass square. At q 2 = follows eq.3 for x(0) if 



the inequality a/A c /i 2 » 1 is assumed. However, at |g 2 
from (73): in (74) there are zeros at the points q 2 = 



m 



v 1 



m 2 (74) strongly differs 
m 2 , but not poles, as it should be 
and as it take place in (73). And also the most important at low Q 2 hadronic term x'ifyQ 2 
is absent in (74). 



VI. SUMMARY 

The g 2 -dependence of topological charge density correlator x{q 2 ) (2) in QCD was consid- 
ered in the domain of low q 2 . For the cases of two and three light quarks the values of x(0), 
obtained earlier [f|- |7]] were rederived basing on the low energy theorems and accounting 
of quasi- Goldstone boson (it, if) contributions. No large N c limit was used and it was no 
appeal to the ^-dependence of QCD Lagrangian (except of the proof of absence of anomalous 
commutator in the sum rule (17)). The only concept, which was used, was the absence of 
Goldstone boson contribution as intermediate state in the singlet (in flavour) axial current 
correlator in the limit of massless quarks. In the three light quark case - the case of real 
QCD - the mixing of 7r and r] is of importance and was widely exploited. The q 2 dependence 
of x(g 2 ) was found as arising from two sources: 

• 1. The contribution of hadronic states (besides 7r and if). This contribution was 
determined from the established in IJ (basing on QCD sum rule approach) connection 
of x'(0) with the part of the proton spin, carried by quarks. 

• 2. The contributions of 7r and r\ intermediate states. These contributions were calcu- 
lated by using low energy theorems only. The final result is presented in eq.73. 



16 



Acknowledgments. 

I am thankful to A.V.Smilga for valuable discussions. I am very indebted to J.Speth for 
the hospitality at the Institut fiir Kernphysik, FZ Jiilich, where this work was finished, to 
A.v.Humboldt Foundation for financial support of this visit and to V.Baru for a help in the 
preparation of the manuscript. This work was supported in part by CRDF grant RP2-132, 
Schweizerischer National Fonds grant 7SUPJ048716 and RFBR grant 97-02-16131. 



17 



REFERENCES 

[1] B. L. Ioffe and A. G. Oganesian, Phys. Rev. D 57, R6590 (1998). 

[2] A. A. Belyaev, A. M. Polyakov, A. S. Schwartz and Yu. S. Tyupkin, Phys. Lett. 59B, 
85 (1975). 

[3] G.'t Hooft, Phys. Rev. Lett. 37, 8 (1976); Phys. Rev. D 14, 3432 (1976). 

[4] R. J. Crewther, Phys. Lett. 70B, 349 (1977). 

[5] G. Veneziano, Nucl. Phys B159, 213 (1979). 

[6] P. Di Vecchia and G. Veneziano, Nucl. Phys. B171, 253 (1980). 

[7] H. Leutwyler and A. Smilga, Phys. Rev. D 46, 5607 (1992). 

[8] B. L. Ioffe, Sov. J. Nucl. Phys. 29, 827 (1979). 

[9] D. J. Gross, S. B. Treiman and F. Wilczek, Phys. Rev. D 19, 2188 (1979). 
[10] M. Gell-Mann, R. J. Oakes and B. Renner, Phys. Rev. 175, 2195 (1968). 
[11] B. L. Ioffe and M. A. Shifman, Phys. Lett. 95B, 99 (1980). 
[12] V. A. Novikov et al, Nucl. Phys. B165, 55 (1980). 



18 



